global dyT Ti T0i dyt1p5 
Ttot = Ti+T0i;
dyT(:,2:end-1) = dyt1p5*(Ttot(:,3:end)-Ttot(:,1:end-2));
dyT(:,1) = dyT(:,end-1);
dyT(:,end) = dyT(:,2);